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Abstract 

Non-relativistic QCD (NRQCD) predicts colour octet contributions to be significant 
not only in many production processes of heavy quarkonia but also in their radiative 
decays. We investigate the photon energy distributions in these processes in the end- 
point region. There the velocity expansion of NRQCD breaks down which requires a 
resummation of an infinite class of colour octet operators to so-called shape functions. 
We model these non-perturbative functions by the emission of a soft gluon cluster 
in the initial state. We found that the spectrum in the endpoint region is poorly 
understood if the values for the colour octet matrix elements are taken as large as 
indicated from NRQCD scaling rules. Therefore the endpoint region should not be 
taken into account for a fit of the strong coupling constant at the scale of the heavy 
quark mass. 



PACS Nos.: ll.10.St, 12.39.Jh, 13.25. Gv 



1 Introduction 



Since the discovery of the J/ ip [||] and the T heavy quarkonia decays are one of the most 
interesting laboratories for investigations within the framework of perturbative QCD. In 
particular these bound states of a heavy quark Q and its antiparticle Q have been examined 
to extract the value of the strong coupling constant a s at the scale of the heavy quark mass 

m Q- 

Early theoretical analyses starting from the calculation of the total rates in leptonic and 
inclusive hadronic decays |J were done in the colour singlet model (CSM). It assumes the 
quark-antiquark pair being in the same quantum state n = 2S + 1 j}j ! > on the partonic level as 
the corresponding quarkonium on the hadronic level. In particular the QQ pair has to be 
in a colour singlet state (C = 1) when it annihilates. As a consequence of this requirement 
the underlying partonic process in the radiative decay H — > 7X of a quarkonium H in the 
ground state 3 S\ is the annihilation of the heavy QQ pair into a photon and at least two 
gluons. This process was calculated first in ||]. 

Great theoretical progress in the understanding of bound states of heavy QQ systems 
has been achieved by NRQCD (Non-Relativistic Quantum Chromo-Dynamics) 0. In 
this theory quarkonia decays are factorized into two step processes: the short- distance 
annihilation of a QQ pair with fixed total spin S, orbital angular momentum L, and total 
angular momentum J and its preceding long-distance transition into this state. While 
the partonic subprocess can be calculated perturbatively to definite order in a s the non- 
perturbative subprocess H — > QQ[n}+ soft degrees of freedom is parameterized by NRQCD 
matrix elements. They are the NRQCD counterparts of the wave function at the origin in 
the colour singlet model and give the probability for finding the quark-antiquark pair in the 
quantum state n at the moment of annihilation. In principle the values of these parameters 
are unknown and must be fitted to experimental data or computed on the lattice ]7|. 
Nevertheless NRQCD provides scaling rules which e.g. predict colour octet matrix elements 
being suppressed by powers of the non-relativistic velocity v with respect to the leading 
order colour singlet matrix element. This typical velocity of the heavy (anti) quark inside 
the quarkonium simultaneously serves as expansion parameter of the effective field theory. 
As a result NRQCD describes a decay rate by an infinite sum over matrix elements of four 
fermion operators with Wilson coefficients which on their part are expansions in a s . 

Taking only the leading order in v/c the NRQCD result coincides with the one of the 
CSM. However, subleading terms in the velocity expansion could be still important numer- 
ically: In radiative decays the partonic kernel of a colour octet contribution is enhanced by 
an inverse power of a s {mQ) with respect to the leading order colour singlet mode. While 
the leading term needs two hard gluons in the final state (cf fig. |l](a)) a QQ pair in a colour 
octet state can annihilate into a photon and a single gluon (cf fig. |l|(b)). Thus one may 
also take into account the subleading terms in v/c. 

The photon energy spectrum in the hard subprocess QQ[n] — > 7X of radiative decays 
has been calculated in next-to-leading order perturbative QCD for both the colour singlet 
mode || and the colour octet modes [f|. Another perturbative contribution may become 
important in the upper endpoint region [z = 2E 1 /Mh — > 1) of the spectrum. Due to an 
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Figure 1: Direct contributions to the radiative decay of a QQ pair: on the left (a) one of 
six colour singlet diagrams, on the right (b) one of two diagrams per colour octet mode. 



imperfect cancellation between terms stemming from real and virtual emission of soft gluons 
one could expect potentially large logarithms ln(l — z) to all orders of the perturbation 
theory. A resummation of these logarithms would then give rise to a Sudakov suppression 
~ exp(— a s ln 2 (l — z)). Though an earlier analysis pl| claimed such a Sudakov damping 



factor in the colour singlet mode a more recent work JTTJ] predicts such Sudakov factors in 
the colour octet channels only while the logarithms should cancel order by order a s in the 
colour singlet mode. 

Besides these perturbative contributions several non-perturbative effects have been in- 
vestigated as well. They become important near the phase space boundaries where the 
photon energy fraction z is small or close to 1, respectively. For low values of z there is a 
large fragmentation contribution caused by the collinear emission of a photon from a light 
(anti)quark in the final state. Examples for such processes are diagrammed in fig. |2|. They 
have been investigated in [jl2| for the colour singlet mode and in |J for the colour octet 
channels. 

At the upper endpoint of the spectrum two different sources for non-perturbative effects 
exist. The first one is the phase space effect associated with the hadronization of massless 
gluons into massive final states. This effect usually is considered by a parton shower 
Monte Carlo thereby generating a non-zero invariant mass for the outgoing gluon(s) [p~3j . 
Another method based on the introduction of an effective gluon mass could obtain 
the appropriate phase space suppression by fitting the values of the effective gluon mass 
to data of radiative J/ip |l5j and T |1| |17| decays independently |18| . 

In this article we concentrate on another non-perturbative effect contributing to the 
upper endpoint of the spectrum. In this region NRQCD operators connected to the center- 
of-mass (cms) movement of the QQ pair inside the quarkonium could become significant 
even though they are subleading in the sense of the naive NRQCD power counting 0, . 
It has been shown in |^0| explicitly that the NRQCD velocity expansion breaks down near 
the endpoint. The reason for this breakdown is the kinematical enhancement of the cms 
operators. In the endpoint region the expansion parameter is v 2 /e rather than v 2 where 
e = I—E^/itlq is a measure for the distance from the endpoint. Thus the velocity expansion 
works fine only for photon energies that are significantly further away from the endpoint 
than AEry ~ vtlqv 2 . However, the range of applicability of NRQCD can be extended to 
higher values for E~ by the resummation of an infinite class of operators into so-called shape 
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Figure 2: Fragmentation contributions to the radiative decay of a QQ pair: colour singlet 
contribution (a), colour octet contribution (b) and (c). The non-perturbative subprocesses 
inside the grey boxes are described by the gluon (a) + (b) and the quark (c) fragmentation 
function respectively. 



functions ||20| . Thereafter the shape function improved spectrum holds up to a resolution 



of AE^ ~ rriQV 2 . 

The shape function formalism yields different shape functions for the different quantum 
states n of the quark- ant iquark pair (QQ[n\). In the leading colour singlet case the shape 
function is nothing else than a number namely the corresponding NRQCD matrix element. 
In the colour octet modes the shape functions are real functions reflecting the phase space 
dependence of the soft gluon radiation. Thus the partonic spectrum of the relevant channels 
n = { 1 S'o 8 ' ) , 3 Pj, 3 S[ 8 ^} which is proportional to 5{E 1 —rag) is smeared out to a quite broad 
peak. 

As mentioned above their contribution may be as large as the one from the leading 
order. If we compare the total rates of the leading colour singlet and a subleading colour 
octet term we get 

47T 

Thus it is worth to model the non-perturbative shape functions to estimate their influence 
on the photon energy distribution. For the construction of the model we keep close to a 
model successfully used for shape functions in quarkonium production |2"Ifl . In this case 



the shape functions which originally were defined in [^| have been modeled by soft gluon 



emission from the final state. Inspired by the success in describing the J ftp production in 
B decays and in the photoproduction channel we take over the physically simple picture 
of radiating off soft gluons from the heavy (anti) quark. 

We will proceed as follows: First we will recapitulate the result of underlying partonic 
process QQ[n] -> 7 + g and QQ[n] -> 7 + g + g for n = 2S + l Lf ] E {tfj , 3 Pj 8) , 3 S[ 8) } and 
n = 3 Si respectively. Moreover we will show to what extent fragmentation contributions 
must be taken into account. Afterwards the construction and application of our shape 
function model in the decay mode is given. Finally we discuss the results of the numerical 
evaluation of the semi-inclusive decay T(1S) — >• 7 + light hadrons. 
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2 The partonic calculation 



Within NRQCD the photon energy spectrum in the semi-inclusive decay H — > 7X of a 
quarkonium is represented by the operator product expansion 

^ = Y,C[n)(H\0[n]\H). (2) 

n 

Here the Wilson coefficient C[n] is calculable perturbatively and gives the differential rate 
dT n / dz in the decay of a QQ pair with quantum numbers n into a photon and light hadrons 
X. Note that in conventional NRQCD the photon energy is normalized on the quark rather 
than on the quarkonium mass: z = E^/uiq. 



2.1 Direct contributions 

Equation (0) includes not only the contributions from the direct production of a photon 
(fig. 0(a)) but also the production via fragmentation (fig. 0(b)). We will come to this point 
later. We first deal with the direct channels. The leading term in the non-relativistic 
expansion is the colour singlet mode displayed in diagram |l|(a). It has been calculated 
perturbatively up to 0(a s ) 0. For the sake of simplicity we restrict ourselves on the tree 
level contribution. It is HI 



^dirr 3Q (i) im _ 32 4^m(m Q )^(m Q ) 
^ [bl 1{Z) ~ 27^ 



2-z z(l - z) 
~z~ + (2 - z) 2 



(3) 



1 — z , , ( 1 — z, 
+ 2— ln(l-i)-2^-^ln(l-5 

The subleading terms 0{v 4 ) in the velocity expansion arise from Feynman diagrams 
like fig. |TJ(b) , where 2S + 1 L^ = l S^\ 3 Pq or 3 P 2 . They are perturbatively enhanced in 
comparison to the colour singlet channel (0(a s ) versus 0(a 2 )). Due to their two body 
kinematics the photon spectrum is fixed to a definite energy value: 

Cf [n](z) = ^-1— i- H[QQ[n] - ig ]{2m Q ) 6(1 - z) . (4) 

The spin and colour averaged squares H[QQ[n] — > r )g\(2mQ) of the amplitudes for the 
three relevant channels are (again we take only leading terms in a s ) |§: 

mnnriQ®! l/o ^ 256ir 2 e 2 Q a em (m Q )a s (m Q ) 

H[QQ[Sy] -»• r yg]{2m Q ) = - , (5a) 

u\nn\3 P (sh 1/0 ^ 7Q8ir 2 e 2 Q a cm (m Q )a s (m Q ) 

H[QQ[ 6 P^ '} -> -fg\{2m Q ) = ^— , (5b) 

H[QQ[ Pi '} -> ig}{2m Q ) = . (5c) 
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2.2 Fragmentation contributions 



Let us turn to the fragmentation contributions now. The corresponding processes are 
associated with diagrams like the ones in fig. |2|. There the photon does not stem directly 
from the annihilation of the heavy quark- ant iquark pair but from the fragmentation of 
a gluon or a light (anti)quark in the final state. Nevertheless we can keep the form of 
equation (0) to describe these contributions, since this subprocess is independent from the 
initial state effects parameterized in the NRQCD matrix elements. The Wilson coefficient 
is then obtained by the folding 



of the fragmentation function D a ^ and the coefficient C^ ir [n] which is the perturbative 
part of the NRQCD decay rate of a QQ pair with quantum numbers n into a particle 
a G {g, q, q} and other light degrees of freedom. The integration variable x indicates the 
energy of the particle a normalized on the heavy quark mass: x = E a /mq. 

Since the fragmentation takes place at a scale far below the heavy quark mass tuq it 
can be factorized from the hard subprocess. This is denoted by the factorization scale /x in 
(^). As usual one may derive a renormalization group equation from the \x independence 
of C 7 rag [n] to shift potentially large logarithms between Cf ir [n] and D a ^. 

Though the fragmentation functions are non-perturbative objects a naive estimate for 
their order of magnitude is obtained from counting coupling constants and collinear sin- 
gularities coming up in a perturbative calculation. For that we look at the subprocesses 
highlighted by grey boxes in fig. |]. The coupling of the photon to a light (anti) quark is 
proportional to a em (mo). For -D 9 _> 7 the leading term comes from the kinematic region 
where the photon and the (anti)quark are collinear. Thus one gets a factor \ti{Q 2 /Qq) 
where Q 2 ~ Wq from the phase space integration with a collinear cut-off parameter Q Q of 

order Aqcd- Accordingly one gets a em (mQ)a s (mQ) ln 2 (<5 2 /Qo) f° r the gluon fragmentation 
function. Here the logarithm appears quadratically because quark, antiquark, and photon 
can become collinear simultaneously in this case. 

The logarithms ln(Q 2 /Q 2 ,) could become so large that they could compensate the per- 
turbative a s suppression and thus confuse the perturbation series dramatically. This is 
seen in a easy way from the running of the strong coupling constant. At leading order the 
renormalization group equation yields 



with p = (33 — 2n/)/3 if rif fermions are active. As long as \i is far above a typical 
hadronic scale /Jo ~ Aqcd one can neglect the constant in the denominator. Hence 



i 



C 7 m %P)= E jfc^[n](x,^)D a ^(z/x,^) 



a=g,q,q ~ 



(6) 




(7) 
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Based on this relation we receive fragmentation contributions with magnitudes comparable 
to the direct ones. This become clear if we rewrite the fragmentation functions in the 
leading logarithmic approximation. Then they have the form ||2T 



2 « em (Q 2 



-/.(0 



(9) 



where / a is a phenomenological function of the fraction £ = E 1 /E a of the photon and 
the parton energy. In case of comparable functions /„ for the different partons a the a s 
in the denominator cancels the additional a s in the hard subprocess of the fragmentation 
contributions. 

Instead of using fitted functions for f a we will take the perturbative result for the 
fragmentation functions D a ^. They are known in next-to-leading order [24| but we will 
restrict ourselves on the leading order for consistency with the investigation of the annihi- 
lation subprocess. The scale dependence of D g ^ and D q ^ is given by the leading order 
DGLAP equations g§ 



d 
d_ 

<9/i 2 



-WO 



Z7T 

1 



(10a) 
(10b) 



where e g is the charge of the light quark measured in units of the elementary charge. The 
Altarelli-Parisi splitting functions are 



5^7 



(0 



i + (i-O s 



wo = 2^ 2 + ( 1 -o 2 



Integrating equation ( 10a|) one obtains the quark fragmentation function 



A,^ 7 (£,/i 2 



2tt 



WO In 



A*g(l"O s 



(12) 



The £ dependence of the starting value -D 9 _> 7 (£, for the evolution in the factorization 
scale /i is mainly determined by a logarithm ln(l/(l — £) 2 ) originating from the phase space 
integration It is already separated in (|12"D. The remaining rest term cannot be calcu- 
lated perturbatively. Therefore it has to be modeled, e.g. with a vector dominance model, 
or fitted to data. The latter was done by the ALEPH collaboration in a measurement 
of the 7 + (1 jet) rates for £ > 0.7 |p8[ . They get the best fit for a constant rest term 
C = -l-ln(Ml/(2^))in 



A?^ 7 (£,^ 2 



2tt 



WO ln 



Mi-0 2 



c 



(13) 
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The non-perturbative scale (jlq extracted from data is then 

^ = 0.14lg;*j GeV C = -13.26l 2 3 ;g. (14) 

The non-perturbative piece of the gluon fragmentation function D g ^ cannot be deter- 
mined experimentally. For the sake of simplicity we will set D g _> y (£, /i ) = 0. The leading 
quadratic logarithmic term reads 

i 

^W£,/* 2 ) = ^ j d ^P 9 ^MD q ^h^ 2 )\^Q . (15) 

Since the experimental investigation of the photon spectrum in radiative quarkonia 
decays is limited to z > 0.4 due to large uncertainties for soft photons caused by 7r° decays 
we need the fragmentation functions for £ > 0.4 only. In this region the contribution of the 
gluon fragmentation function is negligible. However, the quark fragmentation into a photon 
is significant for £ ~ 0.4. Furthermore possibly large contributions in the endpoint region 
we are especially interested in are caused by the logarithmic divergence ln(l/(l — £) 2 ). 

At this stage one comment is in order. Since our approximation for the fragmentation 
functions and in particular relation (|5p holds the better the larger the scale /i is, i.e. the 
larger the heavy quark mass rriQ is, the leading log approximation is inaccurate or even 
not reliable for tuq = m c . Therefore we will concentrate our numerical investigation on 
the T decay as long as we do not restrict ourselves on large values for z. 

Finally we need the perturbative results for the coefficient C^ ir [n] in @. Again we 
distinguish between colour singlet and colour octet contributions. The calculation of the 
colour singlet mode (fig. ||(a) without the fragmentation subprocess) is, except for colour 
factors, the same as the decay rate of ortho-positronium |29[. Similarly the coefficient 
C^ lir [ 3 5'| 1 ' ) ] can be extracted from the photon energy spectrum (ffl). The relative factor 

_ E abc (f d abc )(^ d^) Nl - 4 {Nc=3) 5 

gives the ratio of the corresponding colour traces. Inclusive of the coupling constants the 
triple gluon energy spectrum in the decay QQpS'j 1 ' 1 ] — > ggg results in 



C^S[\x) = B F 2 as{mQ) . CfffflKs) . (17) 
eQOc cm (m Q ) 



This formula contains a combinatorial factor 1/3 = 1/3! : 1/2! which compensates the 
aforementioned factor of three coming from the fact that all three final state gluons can 
fragment into a photon. 

The colour octet coefficients C^"[ 2S+1 L^] are determined by the hard subprocesses in 
diagrams like fig. |2|(b) and fig. 0(c). In both cases the kinematics are trivial. Thus we get 
for Cf r [n] with n G {^ 8) , 3 P (8) , 3 P 2 (8) } from diagram |(b): 

Cf[n]{x) = 2 £-1— -L H[QQ[n] - gg](2m Q ) 6(1 - x) . (18) 
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Up to a factor of two which again is caused by the possibility for both gluons to fragment 
into a photon C^ ir [n](x) matches the corresponding differential decay modes dT n /dx without 
the NRQCD matrix element. The spin and colour averaged square of the amplitudes are 

@: 

HiQQ^] - gg](2m Q ) 

h[QQM 8) } - gg}(2™ Q ) 

H[QQ[ 3 P^} - gg](2m Q ) 

Final states with J = 1 are forbidden by the Landau- Yang theorem |pO|l , i.e. n = 
and n = 3 S'j 8 ' ) do not contribute to QQ[n] — > gg. Instead the latter configuration can decay 
into a light quark- ant iquark pair (cf fig. 0(c)). The corresponding coefficient is 

C?Psf](*) = ± H[QQfs[ 8) ] - qq](2m Q ) 5(1 - x) (20) 

where a G {<?,?} and 

flWgfri"] ^ rf^) ^ 64 ff (m0) . (2D 

Hence we have collected all the results for the partonic decay of a QQ pair that are 
needed for the photon energy distribution in radiative decays of heavy quarkonia. Hence- 
forth we will construct a model for shape functions in quarkonia decays and embed the 
partonic results into it. 

3 The shape function model 

In the following both direct and fragmentation contributions are improved by the applica- 
tion of the shape functions. We model these functions with the following physical picture 
in mind. According to the factorization assumption of NRQCD quarkonium decays are 
divided into two subprocesses as illustrated in fig. [| In the first stage the quarkonium 
H radiates off a cluster of gluons with a momentum k = Yli^i ~ 0(rriQV 2 ). The final 
state of this non-perturbative subprocess is a QQ pair in the quantum state n that anni- 
hilates perturbatively in the second stage. Its momentum is given by Pqq = P — £ with 
P 2 = (2mg) 2 , i.e. the non-perturbative momentum i ~ O(Aqcd) measures the off-shellness 
of the quark-antiquark pair. 

It is important to be careful with neglecting non-perturbative momenta in the hard 
subprocess because the light cone component £ + of the QQ off-shellness is responsible for 
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Figure 3: Schematic representation of direct colour octet contributions to the radiative 
quarkonium decay. 
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the shift of the partonic endpoint mq to the physically correct hadronic value 
| 20| . Therefore we will keep the kinematics exactly here. 

Nevertheless we have to model the radiation of the gluons from the initial state. This 
is done by our shape function |^T) 



f 

J n 



Hi 



dk 2 d 3 k 

2tt (2n) 3 2k° 



{2Tif5\p H + k-P + £)<S> n (k; p H , P) 



(22) 



where <3> n is a radiator function parameterizing the emission of a soft gluon cluster with 
total momentum k. In our ansatz 



$ n (k;p H ,P) 



\k\ b -eM-k 2 o/A 2 n }-k 2 eM-k 2 /A 2 n } 



(23) 



the cut-off parameter A n ~ 0(rriQV 2 ) reflect the expectation that the main contribution 
comes from the ultrasoft region where the energy ko and the invariant mass k 2 of the gluon 
cluster are of order ttiqv 2 and (ttiqv 2 ) 2 respectively. The choice 



(8)i 



(8)] 



2 , ftpF* 
= A[ 3 P (8) ] = A 



: b[ 3 S{ 8) ] 

A[ 3 S[ 8) ] 



0, 
cA 



(24a) 
(24b) 



for the constants in (^) is motivated by the fact that the gluon coupling for a Ml mag- 
netic dipole transition from the quarkonium H to QQ[ Sq ) is proportional to the gluon 
momentum while an El or a double El electric dipole transition to QQ[ 3 Pq g] or <5Q[ 3 5'j 8 ' ) ] 
respectively does not have any k dependence. Furthermore the necessity of at least two 
transitions for n = 3 Si suggests the introduction of a factor c = 1.5 to enlarge the average 
radiated energy and invariant mass in this case. Finally we fix a n by the normalization 
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(2tt)< 



OO OO 

fn(£) = -J^ Jdk 2 j dk yjk 2 - P $ n (k;p H , P) = (H\ O n \H) . (25) 



vOfe* 



3.1 Direct contributions 

After we have determined our model ansatz for the shape functions we proceed with the 
implementation of the hard subprocess. First we deal with the direct contributions. As 
mentioned above in leading order only the colour octet modes are interesting in the shape 
function formalism. Thus we start with the expression 

< ir = E / (^)4 2Af g J^^x(^\P-£-p 1 -px)H[QQ[n]^ 1 g}(P,£,p,,p x ) 

dk 2 rl 3 k 

27(2^ (a)V(PH "^ P + ' ) *" ( ' ;;Pi " P) <26) 

where one can easily recognize the shape_function_in the second line of the equation. In 
( p6|) Mh denotes the quarkonium mass, dp^ and dp x the invariant phase space measures 
of the photon and the hard gluon respectively and finally the spin and colour averaged 
square H n of the hard subprocess amplitude is given by (|]). 

Manipulating eq. ( f26l) we start with integration out the four momenta k and px with 
p\ = which determines the light cone components £\ and £ + by the delta functions to 

t\ = {M H - 2m Q ){M H - 2m Q + 2i ) + i + (2£ - £ + ) - k 2 , (27a) 
1 [Am Q {M H - E 7 ) — M 2 H — 2(M H — 2E 1 )£ + k 2 ] . (27b) 



2E^ 

To make use of these relations we also decompose d 4 £ into its light cone components 

dH = Jdcj) jd£ d£ + ^±Q(£l) (28) 
o 

and rewrite the £q integration into a ko integration by £q = k Q — (Mh — 2mQ). Then the 
integration over the azimuthal angular can be performed trivially because the partonic 
process is independent. Analogously the angular dependence in dp ry is integrated out 

trivially and one gets dp 1 = E^dE^Q(E^)/(4:TT 2 ). Finally we evaluate the equations (|27|). 
In simultaneous consideration of the theta function in (p8|) there arise integration bounds 
for ko where the physical ones are 

(Mh - 2E,) 2 + k 2 M 2 H + k 2 

2(M H - 2E 7 ) ~ k °- ~^2Mh~- (29) 
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From this we can read off an upper bound for the k 2 integration 

k 2 < M H (M H -2E y ) (30) 
while the lower one is determined by k 2 > 0. Introducing the abbreviations 

a = (p H - P 7 )+ = M H -2E 1) (3= (p H - p 7 )_ = M H (31) 
we finally end up with 

a/3 (/3 2 +fc 2 )/(2/3) 

rfr fi dir ^ r dk 2 /■ i i „ r „^ r , l 



/•J/W /"II 1 

27 -•-i?[Q0[nH T ?](%W) • ^<&„(*;p ff ). ( 32 ) 



(a 2 +fc 2 )/(2a) 



This is our master equation for the direct colour octet contributions to the photon energy 
spectrum that take into account shape functions effects within our model framework. Note 
that according to our model the partonic rate depends on 



M Q g{k) = y Mjj - 2M H k + k 2 (33) 

rather than on 2rriQ, i.e. the quark mass in the partonic subprocess is effectively larger 
than rtiQ. The colour singlet contributions are obtained by multiplying (|3[) with the colour 
singlet NRQCD matrix element (H\ C?i( 3 Si) \H). Here the heavy quark mass is set equal 
to M H /2. 

3.2 Fragmentation contribution 

The treatment of the fragmentation contributions is done in the following way: First we 
apply our shape function model on the colour octet contribution dT / dE a (QQ[n) — > ad) to 
extend the reliability of the partonic NRQCD calculation up to higher values of the parton 
energy E a in the QQ rest frame. Afterwards we fold the received spectrum dT/dE a with 
the corresponding fragmentation function -D a ^ 7 : 



Afff/2 

< as _ ^ [^ d lt Da ^/E a ). (34) 



dE~ _ J E a dE n 



Note that the upper bound for the parton energy E a in the quarkonium rest frame is given 
by the hadronic value Mh/2 rather than by the heavy quark mass _E™ ax = rriQ which 
defines the endpoint in the partonic calculation. 

In ([34]) the folding of the gluon and the (anti)quark energy spectrum dT% ir /dE a is 
calculated completely analogous to Cj32|) : 

&P (£ 2 +fc 2 )/(2/3) 
n (a 2 +fc 2 )/(2a) 



11 





5 






4.5 






4 






3.5 






3 












2.5 




"a 






Ph 


2 




1 — 1 


1.5 






1 






0.5 


















sum 

fragmentation 
direct octet 
direct singlet 



A = 300 MeV 




z = 2E 1 /M T 



"a 

Ph 




sum 

fragmentation 
direct octet 
direct singlet 




A = 500 MeV 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

z = 2E n /M T 



Figure 4: Photon energy spectrum in the radiative decay T(1S) 
for two different values of the shape function model parameter: A 
A = 500 MeV (right). 
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where the sum again runs only over the colour octet modes, a and (3 are obtained from 
(PT|) by the substitution — > E a . As described in subsection ^72] the partonic subprocess 

QQ [n] — > g g contributes only for n = ^Sq 8 ^ , n = 3 P > an d n = ^ while n = 3 Si needs a 
light quark-antiquark pair in the final state. The sum in (|34]) runs not only over the gluon 
but also over different (anti)quark flavours q = {u, d, s(, c)} where we additionally assume 
D q ^ ( (0 = D^ ( (0- 



4 Results 

As mentioned above we restrict our numerical analysis on the bottomonium system since 
the fragmentation contributions in the charmonium sector are not very trustworthy within 
our approach. For the photon spectrum in the radiative decay of a T(1S) we choose the 
following set of parameters: Mj- = 9.46 GeV, a s (fi 2 ) = 0.190, a em (fi 2 ) = 1/132 where the 
factorization scale [i is fixed at the b quark mass m;, = 4.8 GeV. Unfortunately the values 
of the NRQCD matrix elements are unknown. While the production matrix elements of 
the bottomonium sector have been fitted in a recent analysis by Braaten, Fleming and 
Leibovich we have to make do with the NRQCD scaling rules 

(T| O s (%) IT) ~ 0HjM|°) ^ (T| 8 ( 3 Sx) IT) ~ v 4 (T| 0^%) IT) (35) 

to estimate at least the order of magnitude of the decay matrix elements needed for the 
normalization, i.e. for the relative weights, of the different contributing channels. For our 
numerical evaluation we take for each colour octet matrix element 2.20 ■ 10~ 2 GeV 2 and 
(T| Oi( 3 Si) |T) = 3.43 GeV 2 , i.e. we set v 2 = 0.08 for bottomonia. Finally we normalize 
the total rate to one. In case of a comparison with data this could be changed easily. 

The result for two different values of our shape function model parameter A ~ m^v 2 
is show in fig. [|. One recognizes that fragmentation contributions are non-negligible for 
small photon energies only. In the region 0.1 < z < 0.4 they are mainly dominated by the 
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Figure 5: Comparison of the theoretical spectrum (solid line) for A = 300 MeV (left) and 
A = 500 MeV (right) with data of CLEO |I7|| . The dashed line shows the direct, the dotted 
line the fragmentation contribution. The spectrum predicted by the colour singlet model 
(direct + fragmentation, without any hadronization model) is indicated by the dash-dotted 
curvature. 



3 S[ channel. Though this channel is connected with the fragmentation function Dg_». 7 (£) 
which diverges for £ — > 1 its numerical contribution to the spectrum for large values of z is 
negligible. This matches to our expectation that fragmentation processes prefer to transfer 
small energy fractions from the gluon (quark) to the photon. 

The upper endpoint region z > 0.8 is dominated by direct colour octet contributions (at 
least in leading order a s ). While the partonic result is proportional to 5(1 — z) the soft gluon 
radiation in the initial state smears out the delta peak and also shifts its maximum to Zqq = 
M^-/Mr < 1. The higher the model parameter A is the smaller is Zqq and therewith the 
effective heavy quark mass and the broader is the width of the peak. Comparing the 
integrated rates of the direct colour singlet mode 

35.(1) 

and the direct colour octet ones 
we realize that their contributions are almost equal. This can be explained by the fact 
that the suppression factor v A = 6.4 • 10~ 3 of the colour octet modes is canceled by an 
additional factor a s (rrib)/(4:n) = 1.5 • 10~ 2 in the colour singlet mode (cf equation ([[])). 
Furthermore one have to consider that there are contributions from many colour octet 
modes ( 2S+1 Lj = 1 Sq , 3 S[ , and 3 P 2 ) but only from one colour singlet mode ( 3 S'{ 1 ' ) ). 

Let us finally concentrate on the region of middle high photon energies. From the 
theoretical point of view the part between 0.4 < z < 0.75 is the cleanest one of the 
spectrum. Here the colour singlet contribution which is assumed to be dominating over 
the whole photon energy range in the colour singlet model can be extracted within shape 
function improved NRQCD without any pollution from other channels. This is still true 
after including the next-to-leading order contributions: While the colour singlet corrections 
steepen the slope the colour octet terms are negligible for 0.4 < z < 0.75 0. 

Since the extraction of a s (m{>) needs an extrapolation of the measured spectrum towards 
z = 0, a complete theoretical understanding of the part of the spectrum the fit is based 



13 



on is indispensable to reach an accurate value for the strong coupling constant. In our 
opinion this is not possible for high photon energies. To illustrate the problems in the 
upper endpoint region we smear out our theoretical result with the energy resolution of 
the CLEO detector 

f = + 1.9-0.125 (36) 

and fit it to their most recent data [lTj for 0.4 < z < 0.7. The result is shown in fig. [5]. While 
the spectrum is described satisfactorily within the fragmentation uncertainties for small and 
middle high values of z the discrepancy between theory and experiment is overwhelming 
for z > 0.75. As shown by the CLEO collaboration |T7| the CSM result combined with 
the fragmentation contributions according to [12] can be brought into agreement with data 
using a hadronization model by Field [T3[] even though the consequential value of a s (irib) 
is slightly to small compared to the one measured on the Zq resonance. 

Although the colour octet contributions in fig. |5| have not been suppressed by a hadroni- 
zation model yet, they seem to be in strong contradiction to the experimental observation. 
The simplest explanation for this deviation would be an extreme smallness of the colour 
octet NRQCD decay matrix elements even smaller than the v A suppression still acknowl- 
edged by the power counting rules. Some hints for such small colour octet decay matrix 
elements also come from the estimation of the a s corrections ||. Furthermore a recent 
analysis of the corresponding production matrix elements yielded smaller values than 
predicted by velocity scaling rules, too. Although the crossing symmetry between these 
production matrix elements and the ones of the decay holds only in leading order pertur- 
bation theory this could also be interpreted as indication for somehow suppressed (or even 
negative?) values of the NRQCD matrix elements. 

Nevertheless the understanding of the upper endpoint region in the radiative decay 
of the T(1S) is not good enough to conclude convincingly that the colour octet matrix 
elements are extremely small. Without having investigated the Sudakov corrections on the 
colour octet contributions and without a better understanding of the hadronization process 
it seems impossible to give a stringent theoretical prediction for z > 0.75. Furthermore 
the experimental investigation of this spectrum suffers from large systematic problems in 
this kinematical regime, too. 

In summary an extraction of a precise value for the strong coupling constant from the 
radiative T decay seems to be impossible as long the upper endpoint region is included 
in the fit. Unless both theoretical and experimental progress concerning the physics in 
the upper endpoint region is achieved a s (mb) should be fitted from the data between 
0.4 < z < 0.7 only. 
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